library(plotrix) 

x <- runif(10, 1, 10)
y <- runif(10, 1, 20)

graf <- function(x,y)
{
  x11() 
  par(mfrow = c(2,2), bty="l") 
  oldpar <- par(bty="n") 
  boxplot(x, y, col=c("purple", "blue"), xaxt="n", ylab="Valores", main="Boxplot", names=c(substitute(x),substitute(y))) 
  mtext(c(substitute(x), substitute(y)), 1, at=c(1,2)) 
  multhist(list(x,y), border="white", beside=TRUE, col=c("purple", "blue"), xlab="Valores", ylab="Ocorrncias", main="Histograma") 
  qqnorm(x, col="purple", pch=16, ann=FALSE) 
  qqline(x, col="purple", lwd=2) 
  par(new=TRUE) 
  qqnorm(y, col="blue", pch=16, xlab="Quantis tericos da \n distribuio normal", ylab="Quantis da amostra", yaxt="n", main="Grfico quantil-quantil")
  qqline(y, col="blue", lwd=2)
  sum.x=summary(x)
  sum.y=summary(y)
  cor.teste=cor(x,y) 
  final <- list(sum.x, sum.y, cor.teste)
  names(final) <- c(substitute(x), substitute(y), "Correlaao")
  return(final)
}

macho <- runif(10, 1, 10)
femea <- runif(10, 1, 20)
graf(macho, femea)


dev.off()





